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£P ■ Abstract 

O \ The results of ATIC-2 baloon experiment (2002-2003) for energy 

spectra of protons, He, C, O, Ne, Mg, Si, Fe, some groups of nuclei, and 
all-particle spectrum in primary cosmic rays are presented in energy 
region 50 GeV-200 TeV. The conclusion is that the spectra of protons 
and helium nuclei are essentially different (the spectrum of protons is 
y—i ■ steeper) and the spectra of protons and heavy nuclei have no-power 

form. 
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<N ; 1 Introduction 

*sO ! 

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment 
designed to measure energy spectra of cosmic rays with individual charge 
resolution from protons to iron in the energy interval approximately from 
50 GeV to 200 TeV. ATIC is comprised of a silicon matrix, a carbon target, 
scintillator hodoscopes, and a fully active BGO calorimeter. The detailed 
description of the spectrometer is given in [TJ, the method of measurement 
£> ! of primary particle charge and the charge resolution achieved in the spec- 

trometer is described in [2]. A matrix of silicon detectors that permits to 
solve the problem of albedo current by the method of high segmentation of 
charge detector was used in the high energy cosmic ray physics for the first 
time [3]. ATIC had two successful balloon flights in Antarctica: from 28 Dec 
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2000 to 13 Jan 2001 (test flight ATIC-1) and from 29 Dec 2002 to 18 Jan 
2003 (science flight ATIC-2). The preliminary results of the processing of 
the ATIC-2 flight data were published in papers [U [5] . 

The calorimeter of the ATIC spectrometer is thin, that is only a part 
of the energy of primary particle is detected by the calorimeter. Therefore 
there exist a problem to reconstruct the spectrum of primary particles from 
the spectrum of energy deposit. The spectrum of energy deposit f(E d ) for 
the particles of some definite sort is connected with the primary spectrum 
$(-Eb) by the Fredholm equation of the first kind: 



where A(Ed, Eq) is the apparatus function which produces for each primary 
energy Eq the distribution of the energies Ed deposited in the calorimeter. 
To obtain the spectrum of primary energies one must solve the equation ([I]) 
for <&(Eq) with known function /(-Ed). This problem belongs to the sort of 
ill-defined reverse problems. Different approximate methods are used to solve 
such problems. The scaling approximation that is used in emulsion calorime- 
ters [6] was used in papers [H E] for transition from the spectra of energy 
deposit to the primary spectra. It is supposed in such approximation that 
the fraction of the energy deposited in the calorimeter (the ratio E^/Eq) and 
the effectiveness of registration do not depend from the primary energy Eq. 
Therefore the primary spectrum could be obtained from the energy deposit 
spectrum simply by energy shift and normalization by the effectiveness of 
registration. 

A simulation of cascades propagation in the spectrometer show that scal- 
ing approximation describes the operation of the device very roughly. The 
method of papers [H [5] should be improved. Two more precise methods 
to solve the problem of reconstruction the primary spectra (CD) are used in 
the present work: direct solution of the Fredholm equation with Tikhonov's 
regularization and the method of differential shifts which is a generaliza- 
tion of the known method [6] . Besides, a new more precise calibration of the 
calorimeter and an improved correction of the thermal sensitivity of calorime- 
ter measurement tract were used in present work. Instead of using of mean 
flight temperature of the device in temperature correction method as in pa- 
pers [HE], instantaneous values of temperature were accounted and used in 
new improved method of correction of thermal sensitivity. 
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2 Reconstruction of primary spectra 



The experimental differential spectra of energy deposited in the calorimeter 
are recorded as a number of counts in logarithmic-equidistant bins that divide 
the whole region of the investigated energies onto n parts. If the primary 
spectra is written in the same form then the integral equation ([T]) in discrete 
form may be written ClS db system of equation 



where the notations are obvious. The search of the primary spectrum 
{Nj,j = 1, . . . ,n} based on the direct solution of discrete Fredholm equa- 
tion ([2]) with using of regularization leads to minimization of some functional 
and using the Monte Carlo method for calculation of statistical errors. The 
realization of the method for processing of the ATIC data was described in 
details in [7j. 

It was shown that the method of solution of the reverse problem with 
regularization for protons and helium nuclei leads to stable result. This 
means that if one processes a spectrum of energy by parts then one obtains 
the respective parts of the primary spectrum that are in the accordance with 
each other. But it turned out that it was not the case for nuclei Z >= 6. So 
the method of differential shifts that generalizes the known method [U] based 
on exact scaling was developed to process the spectra for nuclei Z >= 6. 
The energy scaling factor is calculated separately for each value of Ed in the 
differential shifts method. For the deposited energy value that is exactly 
equal to the logarithmic mean of 2-th bin of the spectrum (Ei) the primary 
energy is reconstructed by the formula 



where K(Ej) is an initial approximation for the primary spectrum. For 
intermediate values of Ed the primary energy is estimated by interpolation 
of the equation ([3]). As an initial approximation a power-form spectra with 
the integral index 7 = 1.6 for momentum were used. 

It was checked that the method of reconstruction of primary spectra by 
solution of reverse problem with regularization and the method of differen- 
tial shifts produce the same (up to statistics) result in the common for both 
methods energy region of the primary spectra. But for all, the method of 
differential shifts is stable for all species, though the lower threshold for re- 
constructed spectra is approximately 200 GeV instead 50 GeV in the method 
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Mi = Y t aaN j , j = l,2,...n, 



(2) 




Sj=i ajjE j K(E j ) 



(3) 
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of regularization. Another advantage of the method of differential shifts upon 
direct solution of reverse problem is that the first is usable in the region of 
data with low statistics where the direct method is already unusable. But it is 
needed for the differential shifts method that the spectrum was not different 
from the power one too much. It is the case for nuclei Z > 6. 

The absolute fluxes of primary particles are determined by normalization 
of the reconstructed primary spectra by geometrical factor of the device, live 
time of the exposition, correction on the residual atmosphere, and correction 
of effectiveness of algorithm of trajectory reconstruction. 

3 Results 

The obtained primary spectra of protons and helium nuclei (here and every- 
where in this paper in the terms of energy per particle) and the ratio of the 
fluxes of protons and helium nuclei against the primary energy are shown in 
fig. [TJ The ATIC-2 results are compared with the results of experiments AMS 
Pi, CAPRICE-98 [ID], BESS-TeV [TT|. The spectra of protons and helium 
nuclei have definitely different slopes and are in the accordance with the data 
of magnetic spectrometers [8] [11] in the region of low energies (where the 
data of magnetic spectrometers are defined). 

The primary spectra of some even nuclei Z > 6 are shown in fig. [21 The 
data of ATIC-2 are compared with the results of HEAO-3-C2 [12J, CRN [T3], 
TRACER (preliminary results) [H] . One can note a good agreement except 
a little systematic exceeding of HEAO data over ATIC-2 data in the energy 
region 200-700 GeV. 

The spectra of different nuclei groups and the spectrum of all particles 
from ATIC-2 data are shown in fig. [3J One can notes a definite flating 
of spectra for nuclei heavier than helium in the region of primary energies 
Eq > 10 TeV. It especially clear in the summary spectrum of nuclei Z > 3 
(fig. [3]) The shape of the spectra of all particles is close to the power low, 
but in the light of the results of ATIC-2 experiment this is seen as a random 
event since the spectra of different nuclei groups (including the spectra of 
protons and helium, fig. [T]) are rather far from a power law. 

The analysis of ATIC-2 data have been continuing. 
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Figure 1: The spectra of protons and helium nuclei (at the left) and the ratio 
of the pronon flux the helium flux (at the right) 
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Figure 2: The spectra of even nuclei: C, O, Ne, Mg, Si, Fe. 
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Figure 3: The spectra of different nuclei groups and the spectrum of all 
particles. The upper threshold in the spectrum of all nuclei Z > 3 is higher 
than in the spectrum of all particles because the spectra of protons and 
helium (see. fig. [T]) can be reconstructed by the method of regularization 
only up to 80 TeV. 
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